Detecting Structure of Complex Network by Quantum Bosonic Dynamics 
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We introduce a non-interacting boson model to investigate topological structure of complex net- 
works in the present paper. By exactly solving this model, we show that it provides a powerful 
analytical tool in uncovering the important properties of real-world networks. We find that the 
ground state degeneracy of this model is equal to the number of connected components in the net- 
work and the square of coefficients in the expansion of ground state gives the averaged time for 
a random walker spending at each node in the infinite time limit. Furthermore, the first excited 
state appears always on its largest connected component. To show usefulness of this approach in 
practice, we carry on also numerical simulations on some concrete complex networks. Our results 
are completely consistent with the previous conclusions derived by graph theory methods. 
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Currently, the study of complex networks has become 
an important research field in physics, biology, sociol- 
ogy, information technology and other branch of sci- 
ences p], H, Hj]. A characteristic feature of these sys- 
tems is the stochastic diffusion of some discrete ob- 
jects on them. Intuitively, such process is referred to 
as granular flows. Perhaps the most well-known ex- 
ample is the flow of information on Internet. In this 
case, messages are encapsulated into discrete date pack- 
ets to be sent from one computer to another. It has 
been observed that, despite randomness of each single 
packet hopping, the motion of massive packets, governed 
by network protocols may reach certain nonequilibrium 
steady states, such as self-similarity of the Ethernet traf- 
fic [J. An important issue arose is whether some gen- 
eral principles can be found to describe dynamics of 
granular flows in a given complex network [5(. Various 
techniques have been developed to address this problem 
0, SJSLM EL HI HI El HE ES S3 ■ For instance, in 
Refs. [lot Il7j. both itinerant fermion and boson models 
were introduced to study the distribution of information 
packets on Internet. 

On the other hand, since granular flow on a complex 
network is completely determined by its topology and 
protocols, one would naturally like to ask whether knowl- 
edge on dynamics of these particles can be conversely 
used to detect structure of the network, such as its con- 
nectedness as well as the number of its components. Pre- 
viously, this problem has been studied in framework of 
percolation theory fill epidemiological processes [lj|, 
and network search [2fJ. For example, the authors of 
Ref. [HI has shown that spectrum of the branching ma- 
trix can be used to uncover properties of the largest con- 
nected component of a network. However, a unified ap- 



proach is still lacking. 

In this paper, we shall introduce an itinerant bosonic 
model, whose hopping amplitude between any pair of 
nodes is determined by the local topology of networks. 
As we show, information on their global properties can 
be derived from dynamics of this quantum system. In 
particular, we find that the ground state of this model 
gives the static distribution of particles and its degener- 
acy accounts for the number of connected components in 
the networks. Furthermore, we are able to identify the 
largest connected component by calculating the energy 
gap between the first excited state and the ground state 
of the system. 

To begin with, let us consider a large but finite com- 
plex network with N nodes and a set of links connecting 
them. Define the the network adjacency matrix A by 
the following rules: If nodes i and j are connected by 
a link, the matrix element Aij is set to be unit; Other- 
wise, Aij — 0. A particle (information packet) can be 
transmitted between two nodes if they are linked. Fur- 
thermore, we assume that the network is undirected. In 
this case, A T = A holds true. With these definitions, the 
degree for a specific node i is given by Ki — An- 

For such a network, one would first like to know 
whether it is connected and which component is the 
largest one, if it is disconnected [22}. To answer these 
questions, we consider a non- interacting boson model, 
whose Hamiltonian is given by 
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where c\ (6j) denotes the boson creation (annihilation) 
operator which creates (annihilates) a spinless boson at 
node i. iy represents the boson hopping amplitude be- 
tween nodes i and j and the summation is over all the 
distinct pairs of nodes. 

To build local topology of the network into this model, 
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we further define f j,- by 



U 




if i 7^ r, 
if i = j. 



(2) 



While the first line in Eq. f2} seems natural, the second 
one demands some explanation. In physics, tn is the local 
chemical potential of particles. It controls the probability 
for a particle to stay at node i: A higher value of ta 
makes such probability lower. On the other hand, for 
a particle hopping randomly in network, large value of 
An^/Ki implies that it has less chance to return to 
node i. Therefore, visibility of the particle at this node 
is decreased. It justifies our choice for the numerator of 
tn. In the meantime, we introduce denominator y/K in 
tu for the purpose of normalization. 

We notice that a similar model has been recently intro- 
duced to studyjocalization of light wave-packet in com- 
plex networks [23[ ■ 

Now, by choosing a natural basis of single-particle 
states \i) = c]|0), we re-write Hamiltonian ([I]) into the 
following matrix 



in 



ART 



H 



-A 2 i 



-A12 

sfKi 



\ 



-A 



JV1 



-A 



N2 



nr. 



\ 



J 



(3) 



where Y^'i stands for Ylijii- We will show that the ground 
state and the first excited state of this model provide us 
with information on the network structure. 

First, we notice that matrix H is semi-positive definite. 
In fact, by making a similar transformation with K = 
diag(y/Ki, y/K 2 , ...^/Kn), we obtain 



H = K^HK 
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To this matrix, we can apply Gershgorin's theorem [2J], 
which tells us that, for a N x N hermitian matrix A, 
each of its eigenvalues must satisfy, at least, one of the 
following inequalities 
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Since a mm — X)'j a «m| holds for each row index m in 
H, we conclude immediately by applying Eq. ([5]) that 
any eigenvalue of the matrix is bounded below by zero. 
Therefore, H as well as H are semi-positive definite. 

On the other hand, a direct calculation reveals that 
vector u = (y/Ki, y/K 2 , •• • , ^/Kn) t is an eigenvector of 
matrix H with eigenvalue A = 0. Therefore, we find that 
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FIG. 1: Number of visitations V vs C in (a) Erdos-Renyi 
networks with N = 1000 nodes and the connecting proba- 
bility p = 0.36 and (b) scale-free networks generated by the 
Barabasi-Albert algorithm, with N = 1000 nodes, and an av- 
erage degree of 6. We assign a random value Wij, which is 
randomly chosen from [1,100], to each ty. We originally set 
M = 5 and N s = 10 10 to ensure the system reach the equilib- 
rium phase. The simulations have at least 10 3 realizations. 



is one of the normalized ground state of Hamiltonian JTJ) . 
According to quantum mechanics, square of coefficient 
Ki j y/^2 n K n equals the probability of finding a parti- 
cle at node i in the ground state. This conclusion coin- 
cides with a well-known result in graph theory that, in 
the static state, the time of a randomly-walking particle 
spending at node i is proportional to Ki [22l | . It indicates 
that the ground state $0 of Hamiltonian ([1]) represents 
actually the static state of the corresponding complex 
network. 

Also, this approach can be easily applied to more gen- 
eral cases, in which nodes are linked with unequal weight 
{wij}. To deal with those cases, we multiply each link Aij 
by Wij and then replace degree Ki by Ci = J2i^i w nAu, 
a quantity which is similar to fitness defined in Ref. [25| , 
in Eq. @ . With these changes, we can show that one of 
the ground state of the Hamiltonian is given by 
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by repeating the above procedure. Therefore, the prob- 
ability for finding a particle at node i should be equal 
to 



Pi 



a. 
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in the infinite time limit. 

Our numerical simulation confirms also these results. 
By displacing randomly M non-interacting particles on 
random scale-free networks with a specific set of linking 
amplitudes, we count the number Uj of visitations to node 
i after each particle hops N steps. As shown in Fig. 1(b), 
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indeed, is strictly proportional to capacity Cj of the 
node. Furthermore, we find that similar results for ER 
networks [Fig. 1 (a)] . 

Up to now, we have not discussed degeneracy of the 
ground state yet. In fact, this quantity gives the num- 
ber of connected components in the network. To make 
this point more clear, let us first consider a completely 
connected network. In this case, matrix H in Eq. (|3|) is 
irreducible. Furthermore, its off-diagonal elements are ei- 
ther zero or negative unit. Therefore, Perron- Frobcnius 
theorem in matrix theory applies [2J]. It tells us that 
the lowest eigenvector of the matrix is unique. In other 
words, the ground state <&o is nondegenerate in this case. 

On the other hand, if the network consists of several 
disjoint components, then Hamiltonian fl} has a nonde- 
generate ground state on each of them. Moreover, as 
shown above, the lowest eigenvalue of the Hamiltonian 
is always zero. Consequently, degeneracy of the ground- 
state energy E = is equal to iV c , the number of con- 
nected components in the network. 

To show usefulness of this result, we investigate numer- 
ically some widely studied complex network, the so-called 
Poisson random graph. In this model, the probability for 
degree of node i taking on a specific value k is given 
by 

where c > is the distribution mean. In particular, we 
consider an ensemble of graphs. Each of them has n 
nodes and m edges and appears with equal probability. 
In the following, we denote this ensemble by G n ^ m . 

Our data are shown in Fig. 2. We find that, for 
m < n/2, degeneracy of the ground-state is roughly equal 
to D = n — m. However, when z = 2m/ 'n > 1, D de- 
creases non-linearly as m further increases. This change 
of behavior around z = 1 can be clearly seen by taking 
the first order derivative of the fitting curves. It sug- 
gests that a dynamic phase transition may occur at point 
2 = 1. Our findings are consistent with the previous 
graph theoretic results [22j |. 

In addition, we observe also that the first excited state 
of Hamiltonian (TT]) is generally nondegenerate in our sim- 
ulations. Moreover, its wave function is actually located 
on the the largest connected component of the network. 
In other words, the coefficients bi's in the expansion of 
the first excited wave function 

<t> 1 = J2bi\i) (10) 

i 

are only nonzero for those nodes which belong to the 
same connected subset of the largest size. Our data are 
shown in Fig. 3. From it, one can easily see that, as m 
varies, the number of nonzero expansion coefficients in 
$! is always equal to the size of the largest connected 
component, which is figured out by the standard depth- 
first search traversing algorithm. Moreover, there always 
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FIG. 2: The number of components in random graphs with 
n vertices and m edges. Each point is an average over 10 3 
networks of 200 vertices each. 



exists a single largest connected component in each real- 
ization of network. 

From quantum mechanical point of view, this obser- 
vation is not difficult to understand. Let us think of 
each connected component of the network as a three- 
dimensional infinitely deep potential well. Then, as is 
well known in quantum mechanics, the energy gap be- 
tween the first excited state and the ground state of 
Hamiltonian (Q]) in each well is proportional to V~ 2 / 3 , 
where V is volume of the well [26|. In other words, the 
gap decreases as the number of nodes in the connected 
component becomes larger. Consequently, one will ex- 
pect that the smallest energy gap is reached in the largest 
well (connected component). On the other hand, we have 
shown above that all the ground state energies on the 
connected components of network are fixed at Eq = 0. 
Therefore, the global first excited state must appear in 
the largest connected component. 

Furthermore, we conjecture that the higher excited en- 
ergy spectrum of this model can be used to reveal other 
characters of networks, such as modularity. We shall ad- 
dress this issue in future investigation. 

In summary, we introduce a non-interacting boson 
model on complex networks to investigate their struc- 
tures in the present paper. Our motivation is based on 
the following observation: It is the structure of a spe- 
cific network that determines the flow of particles in it. 
Therefore, our knowledge on particle dynamics should 
be conversely useful to uncover topology of the network. 
Indeed, we show that the ground state energy of this 
model is degenerate if the network consists of N c disjoint 
connected components and its degeneracy coincides with 
N c . Moreover, the square of each coefficient in the ex- 
pansion of ground state gives the correct averaged time 
for a random walker spending at a specific node i in the 
infinite time limit. Finally, we show also that the first 
excited state of the model is always supported on the 
largest connected component of network. Therefore, it 
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gives us a practical way to detect topological structures 
of complex networks. 
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FIG. 3: The largest component size Si and the number of 
<E>i's nonzero expansion coefficients S2 in ER random networks 
with n — 200 vertices and m edges. Red points indicate the 
number of nonzero expansion coefficients in $1, while blue 
points are the results of the depth-first search traversing al- 
gorithm on the same network. Each point is an average over 
I0 3 networks of 200 vertices each. 
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